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THE INFLUENCE OF THERMAL DISTORTION ON THE PERFORMANCE 


OF GRAVITY-STABILIZED SATELLITES 
By Gerd Kanning 
Ames Research Center 


SUMMARY 


An analysis is made of the effect of thermal distortion of booms on the 
attitude of a gravity-stabilized satellite. The attitude computed when only 
the effects of thermal distortion on the solar pressure torques were consid- 
ered is compared with the attitude computed when the interdependent effects 
of changing geometry and mass distribution were considered as well. 

The calculations were made for several satellites having either symmetri- 
cal or asymmetrical boom arrangements. Calculations for all arrangements 
were made for low orbits inclined 45° to the sunline. In addition, calcula- 
tions were made for one symmetric configuration for synchronous altitude. 

The results indicated that including the effects of changing geometry and 
mass distribution is not important for symmetrical arrangements, but is 
important for asymmetrical boom arrangements. Furthermore, for asymmetrical 
boom arrangements, it was found that the thermal distortion can provide some 
compensation for the unbalanced solar pressure torques attributable to asym- 
metry. Thus, the best pointing accuracy might be achieved with booms designed 
to have a specific thermal distortion. 


INTRODUCTION 


For gravity stabilization of a satellite large differences in the 
principal moments of inertia are required to achieve satisfactory restoring 
torque levels. These differences are usually obtained by erecting one or 
more booms from the main satellite body, or "stabilized package." The booms 
are also a source of disturbance arising from solar radiation pressure forces. 
Careful design of the satellite is required to minimize attitude errors 
resulting from this disturbance. 

The obvious way of minimizing the attitude errors is to minimize the 
unbalanced solar torques. This can be accomplished, in principle, by (a) 
erecting the booms symmetrically and (b) constructing the booms so that they 
do not bend when heated by solar radiation. Condition (a) ensures that the 
centers of satellite and projected area nominally coincide while condition 
(b) ensures that condition (a) holds for all orientations of the satellite 
relative to the sun. 



Condition (a) imposes a considerable restriction on the permissible 
satellite geometry and, therefore, greatly infuences the overall design. 

The configurations that satisfy condition (a) tend to have a large number of 
booms since they must be arranged in symmetrically disposed pairs, implying 
added complication, lower reliability, and possibly greater stabilization 
system mass than asymmetrical arrangements. If the attitude performance 
requirements of the satellite are not too stringent, the reduced complication 
of an asymmetrical boom arrangement may be more important than the decrease 
in attitude errors gained by minimizing solar torques. This is the situation 
for low altitude satellites for which the gravity gradient is relatively 
large and the attitude penalties for unbalanced solar torques more tolerable. 
It may, therefore, be advantageous from the aspect of overall design not to 
attempt to meet condition (a) . 


It is important to recognize here that failure to meet condition (a) may 
invalidate condition (b) as a means of minimizing attitude errors. As an 
illustration, consider the static equilibrium state of the simple, asymmetric. 


single boom satellite shown in figure 1 
changes the magnitude of the unbalanced 
lite, the effect is small compared with 




Figure 1.- Simple asymmetric two-axis gravity- 
stabilized satellite. 


Although the thermal distortion 
solar torques acting on this satel- 
that resulting from the basic geomet- 
rical asymmetry. Therefore, if the 
mass of the boom is small compared 
with the masses of both the tip and 
the stabilized package, the angle a 
between the vertical and the line con- 
necting the tip mass to the stabilized 
package (fig. 1(a)) is only slightly 
influenced by thermal distortion. But 
thermal distortion (provided it is not 
too great) rotates the stabilized 
package so that the earth-pointing 
axis, fixed in the stabilized package, 
is closer to the vertical (fig. 1(b), 
angle B), and pointing accuracy is 
improved. This example demonstrates 
the purely geometrical benefits of 
thermal distortion. 


These geometric effects on asymetrical satellites are not the only 
benefits of thermal distortion. It also changes the mass distribution of some 
satellites so as to reduce pointing errors. Consider the static equilibrium 
state of the simple, asymmetric, two-boom satellite shown in figure 2, whose 
pointing axis is normally parallel to the principal axis of inertia. With 
the solar radiation in the direction indicated, boom A experiences practically 
no solar pressure, whereas boom B is almost orthogonal to the radiation and 
experiences the maximum solar pressure. Relatively large unbalanced solar 
torques therefore act on the satellite. In equilbrium, the principal axis of 
inertia makes an angle with the local vertical, and for a satellite without 
thermal distortion, the pointing axis makes the same angle. On the other 
hand, thermal distortion rotates the principal axis of inertia so that it is 
no longer parallel to the pointing axis. Since thermal distortion has only 
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Figure 2.- Simple asymmetric three-axis stabilized 
satellite. 

a minor influence on the solar pressure torques, it is clear from figure 2 
that the pointing axes must be rotated toward the local vertical. Pointing 
accuracy is again improved. 

The major objection to these arguments is that, except for certain 
special orbits and satellite configurations, static equilibrium never occurs. 
In general, since the thermal distortion varies with position of the satel- 
lite in the orbit, the mass distribution varies continuously and rigid-body 
dynamics do not apply. An adequate dynamical model must represent not only 
the rates of change of inertia but also the rates of change of the angular 
momentum of the distorting portions of the satellite relative to the rigid 
stabilized package. The simple arguments, based on static equilibrium, never- 
theless, give the intuitive feeling that attitude errors can be reduced by 
permitting a controlled amount of thermal distortion. This possibility merits 
investigation in view of recent advances in boom technology, which allow the 
designer much greater control over thermal bending. 

Flexural motions of booms distinct from the thermal distortion described 
above can also be caused by solar heating. Two types of flexural motion have 
been observed, neither of which is considered in this study. One is the 
result of instability caused by steady heating. This instability can produce 
disturbances far exceeding those resulting from the type of boom motion con- 
sidered in this study (see ref. 2 ) . However, recent improvements in boom 
design should eliminate this type of flexural motion (see ref. 3) . The 
second type of flexural motion is caused by sudden temperature changes as the 
satellite moves into and out of the earth^s shadow (thermal twang). For the 
satellites under consideration, the period of these oscillations will be much 
shorter than the period of the librational motion. Therefore, little coupling 
will exist and the librational and flexural modes can be analyzed separately. 

In this paper, the influence of thermal distortion on one symmetric and 
two asymmetric satellite configurations is examined (fig. .3). All three con- 
figurations are stabilized about all axes and are of the coupled, single dam- 
per type (ref. 1) . In both asymmetrical configurations the fixed booms form 
a symmetrical V fixed to the stabilized package at the apex. In one con- 
figuration the V is bisected by the local vertical and points earthward. In 
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(a) 



(a) Vertically 
symmetric. 


[b) Horizontally (c) Symmetric, 

symmetric . 


Figure 3.- Inertial ly coupled, three-axis gravity-oriented 

satellites . 


the other it is 
bisected by the local 
horizontal plane and 
points approximately in 
the direction of the 
velocity vector. These 
satellites are termed 
the vertically and hori- 
zontally symmetric 
satellites , respectively 
(fig. 3 ). The study is 
directed at determining 
the variation of atti- 
tude errors due to solar 
pressure torque alone, 
over a range of values 
of the radius of curva- 
ture of the thermally 
distorted booms. 


It has been customary in past studies of thermally distorted symmetrical 
satellites to assume that the dynamic effects of variable mass distribution 
are small (ref, 4 ) . The argument for neglecting these terms from the equa- 
tions of motion is based on the nominal symmetry of the configuration and the 
symmetry of the distorted shape. The validity of this assumption is also 
examined in this paper. 


NOMENCLATURE 


B damping constant of angular rate damper; made dimensionless by 

dividing by 

b. unit vectors along the axes of principal moments of inertia of the 

^ main satellite body 


C 


c . 

1 



Fm 


spring constant of mechanical spring connected to the damper boom; 
made dimensionless by dividing by 

unit vectors defining orthogonal reference frame in main satellite 
body: C2 lies along the hinge axis of the damper, lies along 

the damper boom in equilibrium, and C3 lies along V3 

time derivative operator with respect to the body frame 

time derivative operator with respect to the inertial frame 
element of mass of damper boom 
element of mass of main body 
force on particle of mass dm 
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inertia matrix o£ the main body 

I j damper boom moment of inertia about the hinge axis 

Oj unit vectors defining right-hand orthogonal reference frame fixed at 

satellite center of mass: 53 is directed toward the center of 

the earth; Oi is in the direction of motion 

R radius of curvature of a boom to its axis when solar radiation is 

normal 


r^ center of mass location of damper 


■-m 


vector from center of mass of the main body to particle of mass dm 


roll, 
pitch, ^ 
yaw 


attitude angles relating Vj frame to Oj frame; roll is a rota- 
tion about the velocity vector; pitch is a negative rotation about 
the angular velocity vector, and yaw is a rotation about the 
earth-pointing axis 


T total torque on main body 


^d/h 

^e/b 


T 


e/d 


U 


^3 


/ • 

b/i 

/ ■ 

d/i 


torque on the damper boom due to the main body 
external torques on the body 
external torques on the damper 
unit dyadic operator 

unit vectors defining orthogonal reference frame fixed in main 
satellite body; in equilibrium, V:j = o^ 

velocity of particle dm; given by 
damper offset angle; angle between 
orbital angular rate 
angular rate of body with respect to the inertial frame 
angular velocity of damper body with respect to inertial space 


dt 

V2 and C2 axes 


ANALYSIS 


The effects of thermal distortion were analyzed by comparing the calcu- 
lated steady-state response when thermal distortion was assumed to influence 
the solar pressure torques alone, and when it was also assumed to cause vari- 
ations in geometry and the mass distribution. The steady-state motion was 
calculated in each case by a digital integration of the equations of motion 
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over a sufficient time to permit transient motion to become negligible. The 
scalar equations that were integrated were expanded from the vector equations 
of motion given in the appendix A. 

In calculating the variation of the mass distribution caused by thermal 
distortion, the following assumptions were made concerning the bending of the 
booms : 

1. Curvature was assumed to be constant throughout the length of any 

boom. 

2. The curvature for a given solar incidence was identical for all 
booms and was proportional to the absorbed energy per unit length. 

3. The plane defined by a curved boom contains the sunline. 

4. A boom reaches its thermal equilibrium position instantaneously 
following a change in illumination. 

The last assumption eliminates any consideration of flexural modes of 
motion . 

The analysis was limited to a study of several specific satellite con- 
figurations. The configurations were related in that each relied solely upon 
gravity for both stabilization and damping. The principal difference between 
the various configurations was the number or orientation of the booms. 

Details of the geometry, mass distribution, and orbits are given in table I, 
and sketches of each are shown in figure 3. 

It was assumed in calculating the attitude motion that disturbances from 
solar radiation were the only source of attitude error. The calculated atti- 
tude errors therefore are not representative of the total error that might 
exist for an actual satellite. Most of the results are for satellites 
assumed to be in a 1200 km orbit inclined 45° to the sunline. For this 
particular orbit, the satellite passes through the earth’s shadow. The 
effects of thermal distortion on a symmetrical satellite in a synchronous 
equatorial orbit at the time of a solstice were also studied. 


RESULTS AND DISCUSSION 


Asymmetrical Satellites 

The boom arrangements for the asymmetrical satellites were chosen to 
satisfy a given set of moment of inertia ratios with no consideration given 
to minimizing solar disturbances or the response to them. Comparisons of the 
pointing accuracies of the two asymmetrical satellites therefore have little 
meaning since either might be improved considerably, and either might be con- 
siderably different if the inclination of the orbit to the sun were changed. 
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The influence of the variation in mass distribution caused by thermal 
distortion on the attitude motion was quite different for the two configura- 
tions studied as is illustrated in figure 4. For these particular results, 
the booms were representative of the nonperf orated, silver-plated, beryllium- 
copper booms which have a curvature when normal to the sunline of about 457 m. 
For the vertically symmetric satellite, the maximum error in earth pointing 
was underestimated by a factor of about 5, and the error in yaw by a factor 
of about 2 when thermal distortion was assumed to influence only the solar 
pressure torque. In contrast, including all the effects of thermal distor- 
tion was not as important in calculating the attitude motion of the horizon- 
tally symmetric satellite. For this satellite, including all the effects 
changed the phase but not the magnitude of the small pitch and roll errors 
Csee fig. 4(b)), The calculated yaw attitude, however, was sensitive to 
inclusion of all of the effects of distortion. Inclusion of all of the 
effects decreased the amplitude of the calculated yaw motion and indicated a 
change in the yaw bias of about -2° . 




All distortion effects included 

Distortion effect on solar 

pressure torque only ‘ 




(a) Vertical symmetric. 


(b) Horizontally symmetric. 


Figure 4.- The effect of thermal distortion on the performance 
of an asymmetrical satellite; R = 457 m. 

In the Introduction it was argued that a specified amount of thermal 
distortion might reduce the attitude errors. This possibility was investi- 
gated by varying the boom curvature, R, due to thermal distortion, and com- 
paring the root -me an- square (RMS) roll, pitch, and yaw errors. These results 
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Figure 5.- RMS pointing error as a 
function of boom curvature. 


•5r 



All distortion effects included 

Distortion effect on solar 


pressure torque only 



Q. 


-.SL 



Figure 6.- Steady-state performance for a 
low altitude symmetric satellite; 

R = 457 m. 


are suimnarized in figure 5 where it is shown 
that minima do exist for a boom curvature 
other than zero. For the two examples given, 
the largest effect of varying the curvature 
is to reduce the RMS roll angle of the verti- 
cally symmetric satellite from 3° for a typi- 
cal nonperforate^ silver-plated boom 
(1/R = 0.0022 m’ ) to nearly zero for a boom 
with nearly zero curvature. For this partic- 
ular example, the engineering choice would 
probably be a boom designed to have zero 
curvature . 

For the horizontally symmetric satellite, 
the minimum RMS errors occur at a rod curva- 
ture of about 0.0014 m The principal 
effect of changing curvature for this satel- 
ite is a change in the yaw error. Detailed 
examination of the time histories from which 
the results of figure 6 were derived shows 
that the minimum RMS yaw error occurs when 
the steady component of the yaw error is 
zero. It is anticipated that a change in the 
orbit relation to the sunline might alter 
this steady yaw bias and the boom curvature 
for which it is zero. 

It must be emphasized that all of the 
results presented are considered to be iso- 
lated examples. The only conclusion that can 
be drawn is that all of the effects of ther- 
mal distortion must be included when the per- 
formance of gravity stabilized satellites 
with asymmetrical boom arrangements is esti- 
mated. In addition, the possibility of mini- 
mizing the errors through specification of 
boom curvature should be investigated. 


Symmetrical Satellite 

It has been assumed in past performance 
studies of symmetrical gravity stabilized 
satellites that the primary effect of thermal 
distortion was its influence on solar pres- 
sure torques (see ref. 4). The argument 
supporting this assumption is that distortion 
will always result primarily in translation 
of the principal axes relative to the stabi- 
lized package and relative to the orbital 
reference axes. 
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high 




All distortion effects included 


Symmetrical gravity stabilized satellites are considered primarily for 
altitude orbits (such as synchronous orbits) where asymmetry cannot be 

tolerated because of large unbalanced 
solar-pressure torques. For this reason, 
the effects of thermal distortion were 
studied for a synchronous satellite as well 
as for a satellite in a 1200 km orbit. At 
the lower altitude, the attitude errors 
are small (fig. 6). The errors in earth 
pointing and in yaw are underestimated by 
only about 0.1° and 0.5°, respectively, 
when distortion is assumed to influence 
solar pressure torque only. Even though 
the difference in the two estimates is 
small, it represents a large fraction of 
the total error. For a synchronous satel- 
lite, the situation is quite different. 

At this altitude the gravity gradient is 
greatly reduced and the ratio of solar dis- 
turbance torque to restoring torque is 
increased because of restrictions placed on 
the mass that can be allotted to the stabi- 
lization system. A corresponding increase 
occurs in the importance of the solar pres- 
sure torque caused by thermal distortion 
compared to the other effects. This is 
illustrated in figure 7 where the calcu- 
lated motion is shown for the satellite 
with the characteristics listed in table I. 
In this instance, the error attributable to 

^ . 1 - solar pressure torque alone is nearly iden- 

Figure 7.- The effect of thermal distortion on. ^ -.7 -ri 

the performance of a symmetric satellite at Ical tO that calculated when all of the 
synchronous altitude. distortion effects are included. 



CONCLUSIONS 


Thermal distortion of the booms of a gravity stabilized satellite 
influences the attitude of the satellite not only by causing solar pressure 
torques, but by changing the satellite geometry and mass distribution. Such 
changes, while relatively unimportant for a symmetrical satellite, can be as 
important for an asymmetrical satellite as the solar pressure torques. Fur- 
thermore, minimum pointing errors for asymmetrical satellites may occur with 
booms that have a nonzero thermally induced curvature. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, July 22, 1969 
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APPENDIX A 


DEVELOPMENT OF THE VECTOR EQUATIONS OF MOTION 
FOR GRAVITY-STABILIZED SATELLITE 


The equations of motion for the gravity-stabilized satellite configura- 
tions studied are developed here in vector form. Consider the sketch in 
figure 8. The force on particle dm is given by 



dV 


m 


•'m = dm 


( 1 ) 


Then the torque on the satellite will 
be : 


^ S ^m ^ ^m 


( 2 ) 


m 


If equation (1) is substituted into (2) 


Figure 8.- A particle of mass dm moving in 
some arbitrary path. 

The time derivative of ^ ^m^ * 


df 


^ (i-m Vm) - 


or 


dV 


r.„ X 


m 


m dt 


d 

dt 


Substitute equation (4) into equation (3) 


m 


E 

^m ^ 

m 


X V 
''m 

^ ^m 


df 


] 

dt 

\) - 

df 

m 

dt 


dV 
n 

dt 


E 


/_ dV^ 

dm % X dt 


m 


(3) 


dV 


m 


m 


X V, 


m 


(4) 


X Vr, 


C5) 


but 


df 


X v„ = 0 


dt m 


then equation (5) becomes 
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( 6 ) 


m 


X \) 


The velocity of the particle with respect to the inertial frame can 

be expressed by the kinematic relation: 


'df 


V. 


m 


'df 


m 


dt 


\dt 

If equation (7) is substituted into (6) 


“b/i "" 


(7) 




T = > dm ;5-- 
dt . 

m ^ 


r„ X 


m \ dt 


'^^m\l - - \l 

V“b/i ^ ’'m/ 


.^dn, 

m 


dt . 
1 




" (d^r\] 3T7 [f™ “ (“b/i - ^m)] 

■* m ^ 


( 8 ) 


If the first term of equation (8) is expanded 

___ ^ . ,.l2, 

\ ^ dm — 


m 


['. ■ ©J -E'-l'. ■ (©, ['. ■ O. 

^ L -* m ' 


( 9 ) 


Expanding the second term of equation (8) yields 

ae [^n> " (“b/i ^ ^»)] ( 


^m • ^m^ - ^mi-m 


) ■ Vi 


where 0 is the unit dyadic operator. The torque equation can now be 
written as 


I'm • TmU - rmfm 


) ■ Vi 


m 

" (^)b * “b/i ■“ [f" ( 

m 


dr 

m 
dt . 


( 11 ) 


The first part of equation (11) is the inertia dyadic for a system of 
mass particles that can be expressed as d(Ibwb/i)/dt . The second part of 
equation (11) represents the rate of change of the angular momentum of the 
distorting portions of the satellite relative to the rigid stabilized 
package . 

^ The torque T can be separated into the external torques on the body 
Tg/b reaction torque T^y^l^ due to the interaction of the main body 
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and damper boom. For the main body equation (11) can then be written as 
"^e/b “ "^d/b " dtT (^b ' %/i) / "'v 

lily 



Likewise a similar expression can be obtained for the damper boom. 



md 


The vector equations (12) and (13) when expanded into scalar form were inte- 
grated by a digital computer until the steady state was reached. 
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TABLE I.- SATELLITE CHARACTERISTICS 


1 * 


Parameter 

Horizontally 

symmetric 

Vertically 

symmetric 

Symmetric 
low altitude 

Symmetric 

synchronous altitude 

Main Body 

Stabilized package 





Mass, kg 

218 

218 

218 

218 

Diameter, m 

.69 

.69 

.69 

.69 

Booms 





Length, m 

30.5 

30.5 

30.5. 

30.5 

Diameter, ra 

.0127 

.0127 

.0127 

.0127 

Mass per unit length, kg/m 

.0250 

.0250 

.0250 

.0250 

Tip mass , kg 

2.269 

2.269 

.729 

2.187 

Reflectivity 

.92 

.92 

.92 

.92 

Angle to vertical reference, deg 

25 

26 

24 

24 

Principal moments of inertia, kg-m^ 





ii 

4754 

4794 

3705 

9235 


4754 

4794 

3705 

9235 

13 

892 

899 

695 

1715 

Damper Body 
Boom 





Length, m 

15 

15 

30 

30 

Diameter, m 

.0127 

.0127 

.0127 

.0127 

' Mass per unit length, kg/m 

.025 

.025 

.0125 

.0125 

Tip mass, kg 

.562 

.562 

.133 

.544 

Reflectivity 

.92 

.92 

.92 

.92 

Angle to plane of main booms, deg 

238.2 

58.2 

238.2 

238.2 

Moment of inertia about hinge, Ij, 





kg-m^ 

155 

155 

116 

301 

Damper parameters 





Damping, B (dimensionless) 

1.059 

1.059 

1.059 

1.059 

Spring, C (dimensionless) 

5.15 

5.15 

■ 5.15 

5.15 

Total Satellite Mass 

225 

225 

i 225 

231 

Orbit 



; 


Period, hr 

1.89 

1.89 

1.89 

24 

Eccentricity 

0 

0 

0 

0 

Angle between sunline and orbit 





normal , deg 

45 

45 

45 

' 23.5 
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